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First order Fermi acceleration at astrophysical shocks is often invoked as a mechanism for the generation of 
non-thermal particles. This mechanism is especially simple in the approximation that the accelerated particles 
behave like test particles, not affecting the shocked fluid. Many complications enter the calculations when the 
accelerated particles have a backreaction on the fluid, in which case we may enter the non linear regime of shock 
acceleration. In this paper we summarize the main features of a semi-analytical approach to the study of the non 
linearity in shock acceleration, and compare some of the results with previous attempts and with the output of 
numerical simulations. 



1. Introduction 

Shock acceleration has been studied carefully 
and a vast literature exists on the topic, includ- 
ing some recent excellent reviews [|l|,^ . Hence we 
do not try to provide here any extensive introduc- 
tion to the problem, but we rather limit ourselves 
to summarize some of the open issues related to 
the backreaction of the accelerated particles onto 
the shocked fluid. The accelerated particles start 
to affect the fluid when their energy density be- 
comes comparable to the kinetic energy density of 
the fluid. In this regime the test particle approx- 
imation breaks down, and the standard results of 
shock acceleration cannot be recovered. The only 
way to have a complete quantitative picture of 
this problem is to use numerical simulations |^ 
H, but it is useful to have a practical analytical 
tool to understand the physics behind the simula- 
tions and also to take into account the non linear 
effects also when these simulations are not avail- 
able, which is unfortunately the rule rather than 
the exception. 

Numerical simulations show that even when the 
fraction of particles injected from the plasma is 
relatively small, the energy channelled into these 
few particles can be close to the kinetic energy of 
the unshocked fluid, making the test particle ap- 
proach unsuitable. The most visible effect is on 



the spectrum of the accelerated particles, which 
shows a peculiar hardening at the highest ener- 
gies. 

The need to have a theoretical understanding of 
the non-linear effects in particle acceleration fu- 
eled many efforts in finding some effective though 
simplified picture of the problem (see |^ for a dis- 
cussion of these alternative approaches). We will 
compare our findings with those of Ref . . 

In the present paper we summarize the results 
widely discussed in where an approach was 
proposed that provides a very nice fit to the re- 
sults of simulations and is in agreement with pre- 
vious analyical calculations. 

2. Non Linear Shock Acceleration 

The distribution function of the particles accel- 
erated at a planar inflnite shock can be written 
in implicit form as ||^: 
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where we put Up = ui -\- 

qs = : and Rsub = 

sion factor at the shock surface 
velocity upstream (downstream)]. Ninj is the 
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number density of particles injected at the shock, 
parametrized here as Ninj = TjNgas,!, where 
Ngas.i is the gas density at a; = 0^ (upstream). 
Particles are assumed to be injected at the shock 
surface with momentum pmj- Eq. tells us 
that the spectrum of accelerated particles has a 
local slope given by Q{p) = -^i^ - Ti^^^fe- 
The problem of determining the spectrum of ac- 
celerated particles is then solved if the relation 
between Up and p is found. 

The thermodynamic properties of the shocked 
fluid are embedded in the usual conservation 
equations, including now the contribution from 
accelerated particles. The mass and momentum 
conservation equations read: 

Pauo = PpUp, (2) 
paul + PgM ^ Ppul + Pg^p + PcB,p, (3) 

where Pg^ and Pg^i are the gas pressures at 
the point x = +oo and x = Xp respectively, 
and PcR,p is the pressure in accelerated parti- 
cles at the point Xp (we used the symbol CR 
to mean cosmic rays, to be interpreted here in 
a very broad sense). In writing eqs. (^) and 
(^) we implicitly assumed that the average ve- 
locity Up coincides with the fluid velocity at the 
point where the particles with momentum p turn 
around to recross the shock. Our basic assump- 
tion is that the diffusion is p-dependent and that 
therefore particles with larger momenta move far- 
ther away from the shock than lower momen- 
tum particles. At each fixed Xp only particles 
with momentum larger than p are able to affect 
the fluid. Since only particles with momentum 
> p can reach the point x = Xp, we can write 
PcR,p = X Ip"^°^ dpp^v{p)f{p), where v{p) is the 
velocity of particles whose momentum is p, and 
Pmax is the maximum momentum achievable in 
the specific situation under investigation. In re- 
alistic cases, Pmax is determined from geometry or 
from the duration of the shocked phase, or from 
the comparison between the time scales of accel- 
eration and losses. Here we consider it as a pa- 
rameter to be fixed a priori. From eq. (|^) we 
can see that there is a maximum distance, cor- 
responding to the propagation of particles with 
momentum Pmax such that at larger distances 



the fluid is unaffected by the accelerated particles 
and Up = uq. Assuming an adiabatic compression 
of the gas in the upstream region, we can write 

Pg,P = Pg.o {^y^ = Pg,o (^)^', where we used 
the conservation of mass, eq. (||). The gas pres- 
sure far upstream is Pg^ = p^u^/ {'^gM'^), where 
7g is the ratio of specific heats (7g = 5/3 for an 
ideal gas) and Mq is the fiuid Mach number far 
upstream. Note that the adiabaticity condition 
cannot be applied at the shock jump, where the 
adiabaticity condition is clearly violated. 

After some algebra (se for the details), the 
conservation equations imply the following equa- 
tion: 
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where Rtot — uq/u2, U{p) — Up/uo, and we put 
VU = dU/dlnp. Solving this differential equa- 
tion provides U{p) and therefore the spectrum of 
accelerated particles, through eq. (|l|). 

The operative procedure for the calculation of 
the spectrum of accelerated particles is simple: 
we fix the boundary condition at p = Pinj such 
that U(j)inj) = ui/uq for some value of ui (fiuid 
velocity at x ~ 0+). The evolution of C/ as a 
function of p is determined by eq. (||). The 
physical solution must have U{pmax) = 1 be- 
cause at p > Pmax there are no accelerated par- 
ticles to contribute any pressure. There is a 
unique value of ui for which the fluid velocity 
at the prescribed maximum momentum Pmax is 
'^Pmam — "o- Finding this value of ui completely 
solves the problem, since eq. (^ provides U{p) 
and therefore the spectrum of accelerated parti- 
cles. Fig. 1 illustrates the distribution function 
(multiplied by p^) for Mach number at infinity 
Mq — 43, Pinj — 10~^TOC and rj = 10~^ and for 
Pmax = lO'^mc (solid line), Pmax = lO^mc (dotted 
line) and Pmax = lO^mc (dashed line), where m is 
the mass of the accelerated particles. The super- 
imposed symbols show the corresponding results 
for the method in m. In Fig. 2 we plotted the 
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Figure 1. Comparison between the prediction of 
our model (lines) and those of Ref. (sym- 
bols). The three sets of curves are obtained for 
Pmax — lO^mc (solid line and crosses), Pmax — 
lO^TOc (dotted line and stars) and Pmax — lO^mc 
(dashed line and diamonds. 



results of our method for another set of parame- 
ters (indicated in the figure) and compared these 
results with the output of numerical simulations 
reported in It can be easily seen that the 
agreement is impressive. 

3. Conclusions 

We report on a novel semi-analytical approach 
to non linear shock acceleration, which improves 
some previous attempts of other authors. This 
method is in good agreement with the previous 
approaches, and is also in impressive agreement 
with the results of numerical simulations on shock 
acceleration. An extensive sets of predictions of 
this approach and a more complete comparison 
with previous results are presented in The 
most important phenomenological consequence of 
the inclusion of the non linear effects in shock 
acceleration is the hardening of particle spectra, 
which may reflect in a corresponding hardening of 
the spectra of secondary particles (photons, elec- 
trons and neutrinos) generated in the interactions 



Figure 2. Comparison between the predictions 
of our model (thick solid line) and the results of 
simulations (dashed line) and the approximation 
in Ref. § (solid hght fine). 



of the accelerated particles with the environment. 
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